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Abstract 

We present an algorithm for the Single Source Shortest Paths (SSSP) problem in H-minor 
free graphs. For every fixed H, if G is a graph with n vertices having integer edge lengths and 
s is a designated source vertex of G, the algorithm runs in 0{n^^^-^~^ logi) < 0(n^-^^^ logL) 
time, where L is the absolute value of the smallest edge length. The algorithm computes shortest 
paths and the distances from s to all vertices of the graph, or else provides a certificate that G is 
not iJ-minor free. Our result improves an earlier 0{n^'^ logi) time algorithm for this problem, 
which follows from a general SSSP algorithm of Goldberg. 

1 Introduction 

The Single Source Shortest Paths (SSSP) problem is the problem of finding shortest paths and, in 
particular, the distances from a specified source vertex to all vertices of a given directed graph. This 
problem is one of the classical and fundamental problems in computer science and has numerous 
applications. 

Dijkstra [4] gave an almost linear time algorithm for the SSSP problem if all edge lengths are 
nonnegative reals. His algorithm runs in 0{m + nlogn) time where m is the number of edges and 
n is the number of vertices, if one uses the implementation from [8J. The situation becomes much 
more complicated when negative edge lengths are allowed. Bellman [2] and Ford [6J gave an 0{mn) 
time algorithm for the SSSP problem where the edge lengths are arbitrary reals. No superlogarithmic 
improvement over this simple algorithm is known. If the edge lengths are assumed to be integers, the 
fastest known algorithm to date is an 0{my/ri\ogL) of Goldberg [TTj, improving earlier algorithms 
of Gabow [9] and Gabow and Tarjan [10], where L is the absolute value of the smallest edge length 
(if L is assumed to be a constant and m is sufficiently large then there are slightly faster algorithms 
based upon fast matrix multiplication techniques). 

A lot of research has been conducted in cases where the input graph belongs to some important 
family of graphs, as the SSSP problem has numerous applications even when the graph is a grid 
or a plane graph. For sparse graphs with m = 0{n) edges, the above algorithm of Goldberg runs 
in 0(n^'^ log L) time. Unfortunately, even if we assume that the underlying graph has bounded 
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degree, we do not know how to do better. For planar graphs, however, better algorithms are known. 
A strongly polynomial 0(n^/^) time algorithm was first given by Lipton, Rose, and Tarjan |14] . 
based on the seminal result on planar graph separators of Lipton and Tarjan [15]. For integer edge 
lengths, this was later improved by Henzinger et al. [12| who gave an 0(n^/^ log L) time algorithm. 
As will be explained below, one cannot directly apply this algorithm to ff-minor free graphs since 
finding an 0{^/n) separator for such graphs already requires 0{n^'^) time with present methods 
[1]. A significant improvement was made by Fakcharoenphol and Rao [5] who gave an O(nlog^n) 
algorithm for planar graphs with arbitrary real edge lengths. Finally, Klein et al. improved this to 
O(nlog^n) [18]. These two algorithms are based on planarity. 

In this paper we consider a much more general class of graphs, the class of ff-minor free graphs. 
A graph G' is a minor of a graph G if G' can be obtained from a subgraph of G by contracting 
edges. A graph is H-minor free if H is not a minor of G. In this paper, we say that a directed 
graph G is i7-minor free if the underlying undirected version of G is iJ-minor free. The classical 
Kuratowski- Wagner Theorem [13\ [23] states that a graph is planar if and only if it has no nor 
i^3^3 minors. (For three different proofs of the theorem, see |23|.) 

Families of if-minor free graphs, for some fixed graph H, are the cornerstone of the seminal 
theory of graph minors developed over the last 20 years, in a series of more than 20 papers, by 
Robertson and Seymour. These families are, to date, the most studied families of graphs in modern 
graph theory. The graph minor theory of Robertson and Seymour culminated, in [21j, with a proof 
of the profound graph minor theorem, also known as the Wagner's conjecture, that states that in 
every infinite set of finite graphs, there is a graph which is isomorphic to a minor of another. One of 
the consequences of this theorem is that for any surface S (whether orientable or not) there is a finite 
set of graphs F{S), such that a graph can be embedded in S (without crossing edges) if and only 
if it does not contain a graph from F(S) as a minor. (This result actually follows from a restricted 
version of Wagner's conjecture which was already proved in [20J.) For a very recent survey of the 
theory of graph minor see Lovasz |16j . 

Classes of i^-minor free graphs are much more general, however, than the class of planar graphs 
or classes of bounded genus graphs. For example, the class of i^s-free graphs contains all the planar 
graphs, and many other graphs, but there is no bounded genus surface on which all the graphs from 
this family can be embedded. 

The question we try to answer in this paper is the following: Can we obtain an algorithm whose 
running time is significantly faster than the O(n^'^logL) time algorithm of Goldberg, when applied 
to i?-minor free graphs (it is a well-known old result of Mader [T7| that i/-minor free graphs have 
m = 0{n) edges). Our main result yields a positive answer. 

Theorem 1.1 Let H he any fixed graph. For a given H-minor free directed graph G whose edge 
lengths are integers, and a designated source vertex s G V{G), there is an 0{n^^^'^~'^ log L) < 
0(^n^ '^^'^ log L) time algorithm that computes a shortest path and the distance from s to each vertex 
ofG. 
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It should be noted that if G contains an //-minor, the algorithm may still work as designated, but 
is allowed to fail. In case of failure, the algorithm produces a certificate for an //-minor in G. An 
//-model in (the undirected version) of G is a set of disjoint connected subgraphs {X^ : v € V{H)} 
indexed by the vertices of H, such that for every edge uv S E{H), there is an edge xy G E{G) with 
X G Xu and y G X^. Clearly G has an //-minor if and only if G has an //-model. Thus, in case of 
failure, the algorithm produces an //-model of G. 

Our algorithm is based upon the "four steps" algorithm of Henzinger et al. |l2j, but with two 
major modifications, that are required due to the following obstacles. It should be noted that 
|12j works not only for planar graphs, but also for every graph which satisfies an 0(y^)-separator 
theorem (see the next section for an exact definition), assuming that an 0(-yn)-separator can be 
obtained in linear time. Unfortunately, although //-minor free graphs satisfy an 0(-yn)-separator 
theorem, the fastest algorithm for finding such a separator, due to Alon et al. [Ij runs in 0{n^'^) 
time, and hence we cannot use it directly. It has recently been shown by Reed and Wood [H] that an 
0(n^/^)-separator for //-minor free graphs can be found in linear time. Using this result, it is already 
possible to obtain an algorithm whose exponent is smaller than 1.5, but not as good as the exponent 
in Theorem II. 1[ A bottleneck in the algorithm of [12] is the need to solve many all-pairs shortest 
paths (APSP) problems in small, and very dense, pieces of the graph, and each such computation 
requires cubic time in the size of the pieces. We show that by using separators that are slightly larger 
than optimal, we can solve these APSP problems faster, resulting in a significant improvement in 
the overall running time. The main idea is to first create the pieces quickly using large separators 
by using the algorithm of [1^, and within each small piece, use small separators, via the algorithm 
of [1], to enhance the APSP computations. 

We note that recently, Tazari and Miiller-Hannemann [22] obtained a linear 0{n) time algorithm 
for SSSP in //-minor free graphs where edge weights are non-negative reals, extending another linear 
time algorithm of Henzinger et al. [12] for planar graphs with arbitrary non-negative real weights. 
The latter algorithm of Henzinger et al. (unlike their algorithm in the case of negative edge weights) 
works in linear time also for every bounded degree graph that has an 0(n^~'^)-separator theorem. 
Thus, using the Reed- Wood result mentioned above, it directly works for //-minor free graphs with 
bounded degree. It is trivial to transform any planar graph to another planar graph with bounded 
degree, while maintaining the shortest distances. This, however, is far from trivial in the case of 
//-minor free graphs, and Tazari and Miiller-Hannemann cleverly overcome this obstacle. 

The rest of this paper is organized as follows. The next section contains definitions and lemmas 
that are needed for the proof of Theorem 11.11 The algorithm proving Theorem 11.11 is described in 
Section 3. The final section contains some concluding remarks. 

2 Definitions and Lemmas 

A separation of a graph G is a pair {A,B) of vertex sets A,B CI V{G) such that AU B = V{G), 
and there is no edge with one endpoint m A \ B and the other endpoint in B \A. The set ACi B is 
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called a separator of G. Suppose now that w : V{G) — > is a weighting of V{G). For X C V{G) 
define w{X) = J2vex '^i'^) ^^'^ define w{G) = w{V{G)). We say that a graph G with n vertices has 
an (f{n), a) -separator if there is a separation {A,B) with |yl n -B| < /(n), \B) < aw{G) and 
\ j4) < aw{G). We say that a hereditary family of graphs (a family closed under subgraphs) 
satisfies an {f{n), a) -separator theorem if every weighted graph with n vertices belonging to the 
family has an {f{n), a) -separator. 

By the seminal result of Lipton and Tarjan [15], planar graphs satisfy an (0{^/n), 2/3)-separator 
theorem. In fact, they also show how to compute an (O(y^), 2/3)-separator in linear time. Subse- 
quently, Alon, Seymour and Thomas [T] extended the result of Lipton and Tarjan to i/-minor free 
graphs. The running time of their algorithm is 0(n^'^) for every fixed H. 

Clearly, if g{n) > f{n) then having an (/(n), a)-separator implies having a (g(n), a)-separator, 
but maybe the latter can be found more quickly. We thus say that a hereditary family of graphs has 
an {f{n), a, t{n)) -separator algorithm if it satisfies an (/(n), a)-separator theorem and an {f{n),a)- 
separator can be constructed in 0(t(n)) time. We can therefore state the result of [T] as follows. 

Lemma 2.1 For any fixed graph H , the family of H-mirior free graphs has ari {0{y/n), 2/ 2>,0{v}-^))- 
separator algorithm. 

In a recent result, Reed and Wood [19] generalize the result from [Ij in an interesting way. They 
show that a separator for an //-minor free graph can be found more quickly, if we are willing to 
settle for a larger separator. Quantifiably, their result can be stated as follows. 

Lemma 2.2 Let 7 G [0, 1 /2] be fixed and let H be a fixed graph. The family of H -minor free graphs 
has an {O {n^^^ ^ ^ ), 2/ 3, 0{n^^"'))- separator algorithm. Furthermore, if an input graph is not H- 
minor free then an H -model asserting this fact is produced in 0{n^^'^) time. 

Notice that the case 7 = 1/2 of Lemma 12.21 degenerates to Lemma 12.11 

Suppose G is a graph and F C E{G). The region induced by F is the set of vertices incident 
with an edge of F. A partition of E{G) into k parts defines a set of k regions. We say that a vertex 
of some region is a boundary vertex if it belongs to more than one region. Otherwise, the vertex is 
called internal. An (r, s)-division of an n-vertex graph G is a partition of E[G) into 0{n/r) parts, 
so that each region contains at most r vertices and 0{s) boundary vertices. 

Fredrickson [7] showed that for every r, an (r, ■^)-division of an n-vertex planar graph can be 
found in O(nlogn) time by a simple recursive application of the separator algorithm of Lipton and 
Tarjan. His method carries over without change to the more general setting of a hereditary family of 
graphs with an (/(n), a, t(n))-separator algorithm. Thus, for i/-minor free planar graphs we obtain, 
using Lemma 12.21 and [7] : 

Lemma 2.3 Let 7 G (0, 1/2] be fixed and let II be a fixed graph. For any r < n, an H -minor free 
graph with n vertices has an {r,r^'^~'^^/^)-division and such a division can be constructed in 0{n^^"')) 
time. 
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A family of sets Vi,. . . ,Vk is called a delta system if the common intersection of all of them is 
identical to the intersection of any two of them. 

Lemma 2.4 Let G = {V, E) he a directed graph with n vertices and with V = ViU ■ ■ ■ UVk where 
Vi,...,Vk is a delta system of common intersection cardinality t, and k is a constant. Suppose 
that Gi is the subgraph induced by Vi, and suppose that an APSP solution for Gi is given for all 
i = 1, . . . ,k. Then, an APSP solution for G can be computed in 0{n^t) time. 

Proof: We show how to compute the nxn distance matrix of G given the distance matrices of the 
Gi in the claimed running time. (The construction of the nx n predecessor matrix representing the 
shortest paths given the predecessor matrices of the Gi can be computed similarly in the same time.) 

Let T = ni=i,...,jkFi and let Wi = Vi\ T. Thus, V = T U Wi U ■■■ U Wk is a partition of V. Let 
Di denote the given distance matrix of Gi for i = 1, . . . ,k. Thus, Di{u,v) = di{u,v) is the distance 
from n to t; in Gi where u,v £ Vi. Let D denote the distance matrix of G. 

We initially compute D{u,v) in the case where both u,v gT. Define a complete directed graph 
Gt on the vertices of T by setting the edge length of {u, v) to be 

k 

wct{u,v) = min di{u,v). 

1=1 

As k is constant and since the Di are given, the directed graph Gt is constructed in O(i^) time. 

Next, we solve the APSP problem in Gt in 0{t^) time using, say, the Floyd- Warshall algorithm. 
Let Dt be the resulting distance matrix of Gt- We claim that D{u,v) = Dt{u,v) for all u,v e T. 
Indeed, any shortest path from n to in G is constructed of segments of shortest paths, where each 
segment is a shortest path in some Gi from some vertex a; G T to some vertex y e T. Since the 
length of this segment is (at least) the length of the single edge in Gt from x to y, the claim follows. 

We next compute D{u,v) where u £ T and v G V\T. Suppose v G Wi. Either there is a shortest 
path in G from n to f that is entirely contained in Gi, or else there is a shortest path formed by the 
two segments from u to some z £ T and from z to v, where the latter is entirely contained in Gi. 
Notice that the distance of the first segment is already computed as Dt{u, z) and the distance of the 
second segment is Di{z,v). It follows that 

D(u,v) = min Dt(u, z) + Di(z,v). 

We therefore get that all the values D{u, v) where u eT and v eV\T can be computed in 0{nt^) 
time. Similarly, we compute D{v,u) where u eT and v eV\T. 

Finally, we compute the remaining D{u,v) where u,v G V \ T. Either there is a shortest path 
from u to V that contains a vertex of T, or else u, v belong to the same Wi and there is a shortest 
path from u to v that is entirely contained in Gi. Thus, setting D'{u, v) = min^g^ D{u, z) + D(z, v) 
we have that, in case u and v are in distinct Wi then D{u,v) = D'{u,v). In case D(u,v) are in the 
same Wi then D{u, v) = va.m.{Di{u, v) , D'{u, v)}. We therefore get that all the values D{u, v) where 
u,v £V\T can be computed in 0{n'^t) time. ■ 
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3 Proof of the main result 



Throughout this section we assume that H is any fixed graph, G is an n- vertex directed graph with 
integer edge lengths, s S V{G) is a designated source vertex, and —L is the smallest edge length 
appearing in G. We show how to compute the distance from s to each vertex of G in the time stated 
in Theorem II. 1[ The computation of the actual shortest paths (in the form of a predecessor tree) 
will be evident from the description. We also assume that G has no negative length cycles reachable 
from s. The algorithm can be easily modified to detect such a cycle if at least one exists. 

Let l/2>7>0bea fixed parameter to be chosen later, and let r be a function of n to be chosen 
later. We follow the four steps algorithm from [12] and apply the lemmas from the previous section 
in the appropriate places. 

First step: We apply Lemma 12.31 and obtain an (r, r(-^^''')/^)-division of G in 0(n^~^'^) time. We 
obtain a set TZ of 0{n/r) regions, where each R & TZ has \R\ < r vertices, and has boundary B{R) 
with \B{R)\ = 0(r(2-7)/3). 

Second step: The step applies the following procedure to each region R £ TZ. As in [12], the goal 
here is to obtain an auxiliary graph Hfi with the following properties. The vertex set of Hf( is B{R), 
and for each ordered pair of vertices u,v €z B(R) there is an edge {u,v) in whose length is the 
distance from n to v in R. In particular, is a complete directed graph (possibly with some edges 
having infinite length). However, the way we construct Hn is different from that in [12] since we 
must avoid the naive (say, Floyd- Warshall) application of APSP on a graph of size as 
this is too time consuming. 

We will show how to create Hji in 0{r^^^^^'^"'^^) time. Using Lemma 12.11 we find an 0{^/r) 
separator for R that breaks R into three pieces Ti, T2, so that each piece contains a 2/3- 
fraction of R and a 2/3- fraction of B{R). (by "breaks" we mean that after removing X/j from R 
the remaining vertices are partitioned to Ti, T2, and T3). Such a separator can be obtained by first 
finding an 0(\/r)-separator that breaks R into two pieces, each containing at most a 2/3-fraction of 
R, and then finding an 0(y^)-separator of the part that contains more than half of the vertices of 
B{R), so that no more than a 2/3-fraction of the elements of B{R) remain in a part after removing 
this second separator (note that here we use weighted separators) . 

Now define Ri = TiU for i = 1,2,3 and notice that the boundary of Ri is (contained in) 
Xr U {B(R) n Ti), and thus define B{Ri) = Xr U {B(R) n Ti). Consider next each Ri as a new 
(smaller) region, and apply the procedure recursively to obtain auxiliary graphs Hji^ for i = 1,2,3. 
How do we compute Hji, given the recursive computations of the Hr.? Notice that the union of the 
Hji- is a delta system with common intersection Xpt whose cardinality is t = 0{^/r). Furthermore, 
the Hjj- 's are complete directed graphs and the edge lengths of each Hr. constitute an APSP solution 
for the Hfi-. Thus, an APSP solution for the union of the Hji. can be computed, by lemma [23] (with 
n = 0{r^'^^'~<^/'^) in the statement of the lemma) in time 
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Since the union of the Hji.^s contains B{R) (and possibly some other vertices of Xji \ B{R), we 
construct Hji by setting, for all ordered pairs u,v £ B{R), the length of {u,v) to be the distance 
from ti to 1" in the union of the Hji. (this distance is given to us from the APSP solution), and this 
distance is clearly also the distance from u to w in R, as required. Notice that the O(-) notation 
should be replaced with O(-) when summing up the running time over all the recursive calls, as 
this summation adds a logarithmic factor to the running time. Thus, the overall time required to 
create the Hji is 0(r^^/®~^'''/^), as stated. The time to create all the for all regions is, therefore, 

0(„^(5-47)/6)_ 

Third Step: We compute distances from s to all the vertices in U/j^t^ -B(i2). This is done by 
replacing each region R with the complete directed graph Hji. (Notice that an edge {u,v) may have 
multiplicity now, since both u and v may appear together in more than one B{R). Notice also that 
we assume that s belongs to some B{R) (when performing the first step, we have at least one region 
R for which s E -R, so we can artificially add s to B{R), if it is not already there). The number of 
edges in the replaced graph (with multiplicities) is 

The number of vertices in the replaced graph is 0{jA'^^^^^'^). The smallest possible edge length in 
the replaced graph is trivially not smaller than —nL. By applying Goldberg's 0{^/N M log K) single 
source shortest path algorithm for graphs with vertices, M edges, and smallest weight —K on our 
replaced graph, the distances from s to all the vertices in [Jr£TiB{R) can be computed in time 

(5(^3/2^(l-57)/6 iQg^)^ 

Fourth Step: We remain with the need to compute, for each region R, the distance from s to the 
vertices of i? \ B{R). For this purpose, we construct an augmented graph Gr which is obtained 
from R by adding edges, and preserves the distances. That is, if u,v £ R then 6ji{u,v) = 5gji{u,v). 
However, the important feature of Gn is that each distance can be obtained via a path of Gr that 
has only O(logr) edges. Cohen [3] exhibits an efficient construction of such an augmentation, but 
one can also use the construction described in the conference version of [12]. Indeed, in the second 
step, while recursively creating Hr from the Hr-^s, we can also create Gr recursively from Gr-'s be 
defining Gr to be the union of Hr and the Gr.^s. An inductive argument shows that distances of 
R are preserved in Gr, and that each distance in Gr is obtained via a path consisting of O(logr) 
edges (corresponding to the depth of the recursion). Furthermore, the recursion shows that Gr 
has 0(r(^~^''')/^ log r) edges (indeed, recall that Hr is a complete directed graph with 0{r^'^~'^'^^'^) 
vertices) . 

Having constructed Gr, we obtain G'^ by adding a new source vertex sr with an edge from sr to 
each vertex of B{R) whose length is the distance from s as computed in the third step. Computing 
SSSP from sr in using Bellman-Ford requires only O(logr) = O(logn) iterations, and each 
iteration is linear in the number of edges of G'^ which is 0(r(^~^''')/^ logr). Thus, shortest paths 
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in G'^ are computed in 0{r^'^~'^'^'^f^) time for each region, and in 0{nr^^~'^^^/'^) time for ah regions. 
Clearly, the the computed distance from sr to a vertex v £ R \ B{R) equals the distance from s to 
V in G. 

Considering all four steps, the overall running time of the algorithm is 

0(max{n^+^ , nr^^"^^)/^ , n^/V^^"^^)/^ log L , nr(i-27)/3|) ^ 

6(max{n^+^ , nA^-^^^^^ , n^/V^^-^T^/^ log L}). 

For a given fixed 7, the optimal choice for r is n^/(^~^'y^ (in fact, it is n^^^^~^'^\logL)^^^^^"'\ but we 
ignore this negligible improvement in the exponent of log L) which now means that the running time 
is 

13-27 

0(max{n^+'^ , n »+^-< logL}). 
Optimizing with 7 = — 3 < 0.392 the running time of the algorithm is 

0{n'^-^ log L) < 0{n^-^^^ log L) 

as required. ■ 
It is interesting to note that the proof of Theorem 11.11 yields a non-trivial complexity bound also 
when applied to multiple sources. Indeed, the first and second step are not affected by the number 
of sources. The third and fourth step can be applied to each source separately. Thus, if we wish to 
compute distances from a set of sources to each vertex of the graph, the time required is 

0(max{nl+^ , nr(5-47)/6 ^ ^3/2+a^(l-57)/6 ^ ^ „l+a^(l-27)/3|)_ 

It is easy to see that for small values of a this is significantly more beneficial than just performing the 
whole algorithm separately from each source, and is faster than any other presently known method. 

4 Concluding remarks and open problems 

Our algorithm for single source shortest paths in i?-minor free directed graphs uses varying separator 
sizes, and utilizes a tradeoff between the size of and the complexity of finding a separator. It would 
be interesting to find other applications of this method. One such result, for the maximum matching 
problem, appears in [25], although in that result the size of the separator is not varying. 
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